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ABSTRACT 


Aims. We explore different scenarios to explain the chemical difference found in the remarkable giant-giant binary system 
HD 138202 + CD—30 12303. For the first time, we suggest how to distinguish these scenarios by taking advantage of the exten- 
sive convective envelopes of giant stars. 

Methods. We carried out a high-precision determination of stellar parameters and abundances by applying a full line-by-line differ- 
ential analysis on GHOST high-resolution spectra. We used the FUNDPAR program with ATLAS 12 model atmospheres and specific 
opacities calculated for an arbitrary composition through a doubly iterated method. Physical parameters were estimated with the 
isochrones package and evolutionary tracks were calculated via MIST models. 

Results. We found a significant chemical difference between the two stars (A[Fe/H] ~ 0.08 dex), which is largely unexpected consid- 
ering the insensitivity of giant stars to planetary ingestion and diffusion effects. We tested the possibility of engulfment events by using 
several different combinations of stellar mass, ingested mass, metallicity of the engulfed object and different convective envelopes. 
However, the planetary ingestion scenario does not seem to explain the observed differences. For the first time, we distinguished the 
source of chemical differences using a giant-giant binary system. By ruling out other possible scenarios such as planet formation and 
evolutionary effects between the two stars, we suggest that primordial inhomogeneities might explain the observed differences. This 
remarkable result implies that the metallicity differences that were observed in at least some main-sequence binary systems might be 
related to primordial inhomogeneities rather than engulfment events. We also discuss the important implications of finding primordial 
inhomogeneities, which affect chemical tagging and other fields such as planet formation. We strongly encourage the use of giant- 
giant pairs. They are a relevant complement to main-sequence pairs for determining the origin of the observed chemical differences 


in multiple systems. 


Key words. stars: abundances — stars: atmospheres — binaries: general — planetary systems 


1. Introduction 


Most binary systems are thought to have formed from a com- 
mon molecular cloud that shared the same chemical compo- 
sition. This is supported by observations of binary systems in 
star-forming regions (e.g. Reiputh et al. 2007; Vogt et al. 2012; 
King et al. 2012) and by numerical models of binary formation 
(e.g. Kratter 2011; Reiputh & Mikkola 2012). However, a num- 
ber of works studied their chemical composition and found slight 
but noticeable differences between the components of some 
binary systems (e.g. Gratton et al. 2001; Desidera et al. 2004, 
2006). 

It is difficult to explain the origin of the chemical differences, 
and a number of possible explanations emerge. For instance, 
Meléndez et al. (2009) showed that the Sun is deficient in refrac- 
tory elements when compared to solar twins. They suggested 
that the missing refractories might be locked up in terrestrial 


planets. They noted that the solar depletion pattern might be 
explained by removing the combined mass of the terrestrial 
planets from the convective zone alone. Some studies followed 
this idea and tried to explain the differences observed in main- 
sequence binary systems where only one component hosts a 
known planet (e.g. Ramirez et al. 2011; Saffe et al. 2015) or 
a debris disk (Saffe et al. 2016). Alternatively, Booth & Owen 
(2020) suggested that a giant planet can block a mass of dust 
exterior to its orbit at the time of planet formation, preventing the 
dust from accreting onto the star. In this way, these two scenar- 
ios (Meléndez et al. 2009; Booth & Owen 2020) attributed the 
chemical differences observed in binary systems to the planet 
formation process. 

Planetary engulfment events also appear to be a promising 
scenario to explain the observed differences in main-sequence 
binary systems (see, e.g. Gratton et al. 2001; Saffe et al. 2017; 
Oh et al. 2018; Nagar et al. 2020; Spina et al. 2021; Jofré et al. 
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2021; Flores et al. 2024). We note that chemical inhomo- 
geneities detected in the open clusters M67, Hyades, and 
Pleiades were also ascribed to ingestion events (Oh et al. 2017; 
Ness et al. 2018; Spina et al. 2018). In particular, Spina et al. 
(2021) performed a statistical study of 107 main-sequence 
binary systems and suggested that planet engulfment events are 
common, with an occurrence >20—35% in solar-like stars. In 
stark contrast, Behmard et al. (2023) claimed that engulfment 
signatures are rarely detected (a rate closer to ~2.9%) and ruled 
out practically all previous engulfment detections in 10 differ- 
ent binary systems. Instead of ingestion events, Behmard et al. 
(2023) suggested that the differences observed in binary systems 
may be attributed to primordial differences between stellar com- 
ponents. We also note that on more theoretical grounds, chemical 
signatures of planet ingestion seem unlikely (Théado & Vauclair 
2012) because the ingested material is rapidly diluted by ther- 
mohaline mixing. Thus, it is not entirely clear that engulfment 
events can explain the chemical differences observed in binary 
systems. 

Primordial inhomogeneities were also suggested based on a 
study of binary systems by Ramirez et al. (2019) and Liu et al. 
(2021). The two studies found that binaries with a higher 
projected separation d present a higher metallicity difference 
|A[Fe/H]|, and they attributed this correlation to a primordial 
chemical inhomogeneity. However, Ramirez et al. (2019) cau- 
tioned that engulfment events could blur the proposed corre- 
lation, while Liu et al. (2021) noted that atmospheric diffusion 
could also add a possible dependence on the stellar parameters 
for main-sequence stars. We also note that diffusion effects could 
produce chemical differences between stars in open clusters (e.g. 
Liu et al. 2019; Souto et al. 2019; Casamiquela et al. 2020). The 
studies above illustrate that the origin of the chemical differences 
observed in binary systems is strongly debated. In particular, it is 
hard to distinguish possible superposed effects such as primor- 
dial inhomogeneities, planetary ingestion, and diffusion effects 
(e.g. Liu et al. 2016, 2021; Ramirez et al. 2019; Behmard et al. 
2023; Nissen & Gustafsson 2018). 

Most previous works focused on main-sequence binary sys- 
tems. However, we note that giant stars are also widely used 
to perform detailed comparisons in the chemical tagging (e.g. 
Bovy 2016; Ness et al. 2018; Price-Jones & Bovy 2019). Giant 
stars could present evolutionary effects in species such as C 
and N. Some authors preferred to avoid these particular species 
in the chemical tagging (e.g. Bovy 2016) while others pre- 
ferred to include them with a caution (e.g. Ness et al. 2018; 
Price-Jones & Bovy 2019). In any case, it is important to note 
that giant stars are significantly less prone to diffusion effects 
than main-sequence stars because their convection zones are 
deeper and more massive (e.g. Korn et al. 2007; Dotter et al. 
2017, and references therein). Dotter et al. (2017) suggested the 
use of giant stars (e.g., in the APOGEE survey) to minimize dif- 
fusion effects in the chemical tagging. In addition, it is expected 
that signatures of external pollution (including their related con- 
densation temperature T, trends) decrease significantly in giants 
stars compared to main-sequence stars (e.g. Fischer & Valenti 
2005; Pasquini et al. 2007a; Spina et al. 2021). In other words, 
giant stars are thought to be significantly less sensitive than 
main-sequence stars to diffusion and pollution effects. This 
makes giant-giant pairs ideal targets for testing the chemical 
homogeneity assumed by the chemical tagging to determine the 
origin of the observed differences. For example, it would be 
highly desirable to explore the correlation suggested between d 
and |A[Fe/H]| (Ramirez et al. 2019; Liu et al. 2021) using giant- 
giant pairs, which would avoid the suggested blur of the corre- 
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lation due to diffusion and pollution effects. However, very few 
works have studied the composition of giant-giant pairs in detail 
so far (e.g. Torres et al. 2015). 

We present in this work the first detailed analy- 
sis of the chemical composition of the giant-giant binary 
HD 138202 + CD—30 12303 (hereafter, stars A and B), with a 
bound probability >95% from Gaia EDR3 data (El-Badry et al. 
2021) and a projected separation of 73” (~38 575 au). Surpris- 
ingly, we found a significant chemical difference between the 
two stars (A[Fe/H] ~ 0.08 dex), which is largely unexpected con- 
sidering their insensitivity to planetary ingestion and diffusion 
effects. We consider this a remarkable result for several reasons. 
First, we found a difference between the two components in a 
binary system in which we separated the source of the chem- 
ical differences for the first time (and we attributed them to 
primordial inhomogeneities). This would imply that metallic- 
ity differences should in general be taken with caution because 
they might originate in scenarios that are unrelated to engulf- 
ments. Thus, in order to claim ingestion events, it would become 
mandatory to explore additional evidence such as the study of 
T, trends, the lithium content, and the stellar rotation. Second, if 
it is confirmed that the slight chemical differences we detected 
are primordial (which would be more in line with the results 
of Behmard et al. 2023), it could have important implications 
as follows. It would challenge the main assumption of chemi- 
cal tagging, that is, the chemical homogeneity. It would place 
important constraints on formation models of multiple systems 
(e.g. Bate 2019; Guszejnov et al. 2021) and interstellar medium 
(ISM) mixing (e.g. Feng & Krumholz 2014; Armillotta et al. 
2018). Finally, a primordial difference could severely impact 
planet formation, which might help to explain, for example, why 
very similar stars in wide binary systems present different plan- 
etary systems (e.g. Biazzo et al. 2015; Teske et al. 2016). This 
highlights that the study of giant-giant pairs (as in the present 
work) may emerge as an important complement to the study of 
main-sequence pairs. 


2. Observations and data reduction 


Observations of the giant-giant pair were acquired through the 
Gemini High-resolution Optical SpecTrograph (GHOST), which 
is attached to the 8.1 m Gemini South telescope at Cerro Pachón, 
Chile. GHOST is illuminated via 1.2” integral field units that 
provide the input light apertures. The spectral coverage of 
GHOST between 360 and 900 nm is appropriate for deriving 
stellar parameters and chemical abundances using several fea- 
tures. It provides a high resolving power R ~ 50000 in the 
standard resolution mode!. The read mode was set to medium, 
as recommended for relatively bright targets. The observations 
were taken on May 15, 2023, during a GHOST science verifi- 
cation run, using the same spectrograph configuration for both 
stars. The exposure times were 3 x 200s and 3 x 180 s on targets 
A and B, obtaining a final signal-to-noise ratio (S/N) of ~400 
per pixel measured at ~6000 A in the combined spectra for both 
stars. The standard star 18 Sco was also observed with the same 
spectrograph setup, achieving a similar S/N to use as initial ref- 
erence. The spectra were reduced using the GHOST data reduc- 
tion pipeline v1.0.0, which works under DRAGONS?. This is a 
platform for the reduction and processing of astronomical data. 


! https://www.gemini.edu/instrumentation/ghost 
2 https://www.gemini.edu/observing/phase-iii/ 


reducing-data/dragons-data-reduction-software 
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Table 1. Stellar parameters derived for each star. 


(Star—Reference) Tef logg A(Fe/H) Vturb 
[K] [dex] [dex] [km s~!] 
(A-18 Sco) 4952+49 2.51 0.09 0.067 + 0.014 1.82 + 0.06 
(B-18 Sco) 5007+45 2.57+0.08 -0.015+0.012  1.72+0.05 
(A-B) 4952+30 2.52+0.06 0.083 + 0.011 1.82 + 0.04 
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Fig. 1. Differential abundance vs. excitation potential (upper panel) and 
vs. reduced EW (lower panel) for the case (A-B). The filled and empty 
circles correspond to Fe I and Fe Il. 


3. Stellar parameters and abundance analysis 


The fundamental parameters (Ter, log g, [Fe/H], and vj») were 
derived following a similar procedure as in our previous work 
(Saffe et al. 2018, 2019). We measured the equivalent widths 
(EW) of the metallic lines using the IRAF task splot in the 
stellar spectra. The line lists were taken from works of giant 
stars (Jofré et al. 2015; Soto et al. 2021) with updated labora- 
tory data for some lines (Liu et al. 2014a; Meléndez et al. 2014; 
Bedell et al. 2014). We imposed an excitation and ionization bal- 
ance of Fe lines using the differential version of the FUNDPAR 
program (Saffe 2011; Saffe et al. 2018). This code made use of 
ATLAS12 model atmospheres (Kurucz 1993) together with the 
MOOG program (Sneden 1973), using specific opacities calcu- 
lated for an arbitrary composition. The stellar parameters were 
determined using the solar twin 18 Sco as reference, that is, (A— 
18 Sco) and (B-18 Sco), by adopting (5811 K, 4.42 dex, 0.04 dex, 
1.06 km s~') for 18 Sco. These are the highest-precision parame- 
ters currently derived for this solar twin (Liu et al. 2020). We recal- 
culated the parameters of star A, but used B as reference, that is, 
(A-B), as listed in Table 1° and as shown in Appendix A. The 
errors on the stellar parameters were derived following the pro- 
cedure detailed in Saffe et al. (2015), which takes the individual 
and mutual covariance terms of the error propagation into account. 
Figure | shows the abundance versus excitation potential (top 
panel) and the abundance versus the reduced EW (bottom panel) 
for the case (A-B) with an average uncertainty of 0.02 dex. 

The hyperfine structure splitting (HFS) was considered for V I, 
MnI, Col, and Cul using the HFS constants of Kurucz & Bell 
(1995) and performing spectral synthesis for these species. 
We also derived the Lil abundance using spectral synthe- 
sis with the resonance line 6707.80 A, which includes the 
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Fig. 2. Differential abundances (A-B) as a function of T.. The contin- 
uous lines show linear fits to all species (black) and to the refractory 
species (T, >950 K, magenta). The dotted black line corresponds to a 
null metallicity difference between the stars. 


doublet 6707.76 Å, 6707.91 Å and HFS components. We cor- 
rected the Li I abundance for nonlocal thermodynamic equilib- 
rium (NLTE) effects by interpolating in the data of Lind et al. 
(2009). We adopted ALinrre-Lre ~ 0.18 dex for both stars, imply- 
ing that NLTE effects are not significant for the case (A-B). 
The C/N ratio, which is sensitive to evolutionary effects, was 
derived through spectral synthesis. The C abundances were 
obtained from the CI lines between 7111 and 7116 Å, while N 
abundances were obtained by fixing the C content and then vary- 
ing N in the CN band near ~4212 A. 


4. Results and discussion 


The star A is more metal rich by ~0.08 dex than its companion 
(see Table 1). The condensation temperatures of the elements 
were taken from the 50% T, values derived by Lodders (2003) 
for a solar system gas with [Fe/H] =0. The differential abun- 
dances (A-B) are presented in Fig. 2 as a function of Te. We 
discarded from the fits species with evolutionary effects such as 
C and N, and species that could drive the trends, such as Zn, 
Al and Ba. The black dotted line would correspond to a null 
metallicity difference between the stars. We obtained slopes of 
0.80 + 6.68 x 1075 dex K~! and —1.20 + 13.20 x 1075 dex K`! for 
all species and for the refractory species, respectively. The abun- 
dance dispersion in Fig. 2 tends to increase the uncertainties 
in the slopes, which should be taken with caution. The aver- 
age abundances of the volatiles (T. < 950 K) and refractories are 
0.073 + 0.035 dex and 0.081 + 0.010 dex, respectively. Although 
we find no significant T, trend between stars A and B, most species 
differ by ~0.08 dex in the two stars. In the next subsections, we 
explore different scenarios to explain the observed differences. 


4.1. Scenario of the evolutionary state 


If both giant stars underwent first dredge up (FDU) before, 
as is strongly suggested by the low C/N and Li abundances 
(Salaris & Cassisi 2017), evolutionary effects (<0.01 dex) can- 
not explain the observed difference of ~0.08 dex in Fe. Even 
when we assume that only one star passed through the 
FDU phase (which is highly unlikely), evolutionary effects 
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(<~0.04 dex) are not enough to explain the observed difference 
in Fe. The calculations are extensively discussed in Appendix B. 


4.2. Scenario of planet ingestion 


A metallicity variation due to a planet engulfment is diffi- 
cult to detect in stars with a significant convective envelope 
(e.g. Fischer et al. 2004; Fischer & Valenti 2005; Pasquini et al. 
2007a,b; Spina et al. 2015, 2021). For example, Fischer et al. 
(2004) and Fischer & Valenti (2005) suggested that the higher 
convective envelopes of subgiant stars (compared to main- 
sequence stars) should result in a much lower or null increase in 
metallicity due to external pollution. Moreover, Pasquini et al. 
(2007a,b) estimated that the convective envelope of a main- 
sequence 1 Mo star increases about 30—35 times when the star 
evolves to the RGB, changing from ~0.022 Mo to ~0.77 Mo. 
The authors claimed that a possible higher metallicity (con- 
fined to the superficial layers) would easily decrease to the pri- 
mordial values through the deepening convective zone (CZ). In 
particular, they estimated that an excess of 0.25 dex in a solar 
star would become lower than ~1% in a giant star. Similarly, 
for the case of ingestion in main-sequence stars, Spina et al. 
(2015) showed that the mass of the CZ is a critical parameter for 
determining the [Fe/H] enhancement (e.g., their Fig. 2). More 
recently, Spina et al. (2021) showed that planetary ingestion is 
extremely sensitive to the thickness of the external layer (e.g., 
their Fig. 1), and they concluded that giant stars (and low-mass 
main-sequence stars) are not expected to vary their chemical 
composition because their external layers are extremely thick. 
In other words, giant stars are thought to be significantly less 
sensitive than main-sequence stars to engulfment events. 

These arguments disfavor an engulfment event to explain the 
observed metallicity difference between stars A and B. However, 
we estimated the amount of planetary material that star A would 
need to ingest to increase its metallicity by ~0.08 dex compared to 
star B. A first-order approximation is detailed in Appendix B. We 
estimate that star A would need to have ingested between 11.0 and 
150.0 Mjup of planetary material, depending on the adopted con- 
vective envelope mass and metallic content of the ingested planet 
(see Table C. 1). In particular, the ingestion of planets with masses 
near the lower value of 11 Mjup would require extremely metal- 
lic objects, which are difficult to explain with current planet for- 
mation models (see e.g. the discussions in Leconte et al. 2009; 
Thorngren et al. 2016). In this case, a low convective envelope 
mass for star A of only 35% of their total mass would also be 
required. In addition, the ingestion of massive planetary bodies 
should produce a significant increase in the rotational velocity 
of giant stars (Carlberg et al. 2009, 2012; Privitera et al. 2016; 
Stephan et al. 2020). We estimate an increase in the rotational 
velocity of star A between 8.72 and 37.88 kms“! after the planet 
ingestion (see Table C. 1). However, we find no evidence of signifi- 
cantly different projected rotational velocities between stars A and 
B (vsini4.73 + 0.44km s™! and 5.19 + 0.34kms“!, after quadrat- 
ically subtracting Umacro from Hekker & Meléndez 2007). Star B 
presents an even slightly higher v sini than star A, in contrast to 
what is expected after an ingestion by star A (although we caution 
that the axis orientation could still play a role). 

We also estimated the expected increase in the lithium con- 
tent of star A produced by the planet ingestion (see Table C.1). 
After the engulfment of 11.0-150.0 Mjup of planetary mate- 
rial, we estimate that the lithium abundance of star A should 
increase by ~1.50 to 2.10dex. However, the Li abundances 
are A(Li)4 = 0.76 + 0.04 dex and A(Li)g = 1.25 + 0.04 dex (esti- 
mated with spectral synthesis and NLTE corrected following 
Lind et al. 2009). In other words, the planet ingestion predicts 
a strong increase in lithium of star A, while the observations 
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seem to show the opposite, that is to say, a higher lithium content 
in star B. We note that some works suggested that an ingestion 
event could be ruled out if the same component of a binary sys- 
tem presented both a higher Li content and a lower metallicity 
of the pair (Ramirez et al. 2019; Spina et al. 2021), as observed 
in the stars of the present work. The Li abundances of stars A 
and B also appear to be similar to those of other giant stars 
(see Fig. C.1), while the expected Li value after an ingestion is 
considerably higher. The higher Li abundance of star B com- 
pared to star A might at least partially be explained by their 
slightly higher Tep (see Appendix B). It is also important to note 
that the lithium abundances are strongly modified by evolution- 
ary effects such as the FDU and other still uncertain processes. 
Privitera et al. (2016) suggested using a possible increase in the 
rotational velocities rather than the lithium abundances in order 
to detect a planetary ingestion event in giant stars. 

To summarize the ingestion scenario, an engulfment that pro- 
duced the metallicity difference between stars A and B would 
require very high values of ingested material together with 
extremely metallic planets and a very low convective envelope 
of star A. This is highly unlikely. These extreme values should 
also be accompanied by a significant increase in the rotational 
velocity and lithium abundances, which are not detected. The 
iron and lithium abundances also seem to contradict the abun- 
dances expected after an ingestion event. In addition, no planets 
are detected around stars A and B. Therefore, current evidence 
does not support the ingestion scenario. 


4.3. Scenario of sequestered refractories 


Meléndez et al. (2009) found that the Sun is depleted in refrac- 
tory species when compared to solar twins, suggesting that the 
missing refractories might be locked in terrestrial planets and 
in the cores of giant planets. The depletion pattern is explained 
if the combined mass of terrestrial planets is removed from the 
convective zone of the Sun (Meléndez et al. 2009; Chambers 
2010; Kunitomo et al. 2018; Bitsch et al. 2018). When the star 
evolves off the main-sequence, the massive convective zone in 
the red giant phase erases this chemical fingerprint, as suggested 
by Maldonado & Villaver (2016). This implies that the mutual 
chemical difference found between the giant components in this 
work (where no clear T, trend or planets are detected) can hardly 
be attributed to this effect. 

We also note that the massive convective envelopes of 
giant stars could present an additional advantage for testing the 
original homogeneity of the chemical tagging. We mentioned 
that giants stars are less sensitive to diffusion and pollution 
effects than main-sequence stars. In addition, if planet forma- 
tion produces a possible sequestration of refractories in the most 
superficial layers of main-sequence stars (as suggested by e.g. 
Meléndez et al. 2009; Chambers 2010; Kunitomo et al. 2018; 
Bitsch et al. 2018), this fingerprint should be diluted by the mas- 
sive envelopes of giant stars. 


4.4. Scenario of primordial inhomogeneities 


Some works suggested that the slight chemical differences found 
between the components of main-sequence binary systems could 
be attributed to primordial inhomogeneities (Ramirez et al. 
2019; Liu et al. 2021; Behmard et al. 2023). We suggest that the 
chemical differences found between stars A and B are also pos- 
sibly primordial by previously ruling out other scenarios. The 
large separation of the two stars (~38 575 au) also supports this 
idea. However, no 7, trend is detected between stars A and B. In 
this way, the main finding of this work is the detection of clear 
chemical differences but null T, trends in a giant-giant pair. 
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Further relevant implications of this result are listed below. 

(a) In general, binaries with overall abundance differences 
but no clear T, trends are likely due to primordial inhomo- 
geneities (many cases of chemical anomalies in different liter- 
ature works, such as e.g. Spina et al. 2021). 

(b) This poses a crucial challenge to the concept of chem- 
ical tagging. Primordial inhomogeneties very likely stem from 
an inhomogeneous ISM, and therefore, the origin of stars can- 
not be easily tagged. We derived the 3D velocity difference Av 
between stars A and B and obtained 1.13 kms, which would 
support the binarity (i.e., the conatal nature) of this pair. For 
example, Yong et al. (2023) studied a sample of 125 comoving 
pairs, selected by taking Av < 2.0 km s™!, a cutoff chosen to avoid 
nonconatal pairs. 

(c) Primordial inhomogeneities would place important con- 
straints on formation models of multiple systems (Bate 2019; 
Guszejnov et al. 2021) and on hydrodynamic simulations of ISM 
mixing (e.g. Feng & Krumholz 2014; Armillotta et al. 2018). 
For example, Bate (2019) showed that a lower metallicity results 
in an increased fragmentation in cores, filaments, and disks 
because the cooling rate of dense gas is higher. 

(d) A primordial difference might also affect planet forma- 
tion. This might help to explain, for example, why very similar 
stars in wide binary systems can present different planetary sys- 
tems (e.g. Biazzo et al. 2015; Teske et al. 2016). 

We strongly encourage the study of giant-giant pairs. This 
novel approach might help us to distinguish the origin of the 
slight chemical differences observed in multiple systems. 
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Appendix A: Stellar parameters and abundances 


Table A.1. Chemical abundances for the case (A-B). 


Specie [X/H] + Crot T obs O par 
[Li I/H] -0.490 + 0.031 0.010 0.029 
[C I/H] 0.105+0.030 0.007 0.029 
[N I/H] 0.040 + 0.032 0.010 0.029 
[O I/H] 0.048 +0.058 0.010 0.057 
[Na I/H] 0.065 + 0.033 0.021 0.026 
[Mg I/H] 0.077+0.027 0.010 0.025 
[A1 IH] = -0.024+0.022 0.010 0.019 
[Si I/H] 0.115+0.011 0.010 0.005 
[S 1/H] 0.097+0.042 0.010 0.041 
[Ca I/H] 0.041 +0.022 0.012 0.018 
[Sc H/H] 0.105+0.033 0.025 0.022 
[Ti I/H] 0.063+0.013 0.008 0.010 
[V TH] 0.035+ 0.011 0.006 0.009 
[Cr I/H] 0.045+ 0.013 0.008 0.010 
[Mn I/H] 0.070+0.031 0.010 0.029 
[Fe I/H] 0.083 + 0.004 0.003 0.003 
[Fe W/H] 0.083+0.020 0.012 0.016 
[Co I/H] 0.077+0.013 0.007 0.010 
[Ni I/H] 0.097 + 0.008 0.006 0.004 
[Cu I/H] 0.130+0.031 0.010 0.029 
[Zn 1/H] 0.157+0.036 0.010 0.034 
[Y O/H] 0.083 + 0.038 0.010 0.036 
[Ba ll/H] -0.100+0.031 0.010 0.029 
[Ce I/H] 0.124+0.030 0.010 0.029 
We present in Table A.l the differential abundances 


[X/H] + oror obtained for the case (A-B). For each species, we 
show the observational error Cops (estimated as a / y(n -— 1), 
where œ is the standard deviation of the different lines) as well 
as internal errors due to uncertainties in the stellar parameters 
par (by adding quadratically the abundance variation when 
modifying the stellar parameters by their uncertainties). For 
chemical species with only one line, we adopted for o the 
average standard deviation of the other elements. The total error 
ror was obtained by quadratically adding Cops and © par. 

We also took the opportunity and estimated Teff of stars 
A and B by using different photometric calibrations in order 
to compare them with high-precision spectroscopic results. 
Ramirez & Meléndez (2005) provided the T.¢¢ calibration as 
a function of (B-V) through a polynomial fit. We derived an 
extinction of A, ~ 0.57mag from Gaia Collaboration (2018) 
and A, ~ 0.73mag from the maps of Schlegel et al. (1998). 
By using A, from Gaia, we obtain T.¢¢(A)=49444126K 
and T,¢¢(B) =4991+128 K, while with A, from Schlegel et al. 
(1998), we obtain slighly higher values, T, (A) =5054+132 K 
and T,¢(B) =5104+140 K. We note that the spectroscopic Tess 
agrees well compared to the photometric T.¢r with A, taken 
from Gaia. This indicates that the extinction values from Gaia 
might be more appropriate. By using this photometric estimation 
of T. +, the difference between the two stars would be T,¢ (B) — 
Ter¢(A) ~ 47 + 179 K, which agrees well with the spectroscopic 
difference of 55 + 66 K. However, the photometric estimation of 
the stellar parameters should be taken with caution due to the 
significant extinction in the direction of the stars. Stars A and B 
have a distance of 523+6 pc and 536+6 pc, respectively, based 
on Gaia EDR3 parallaxes. These distant stars could lead to dif- 
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ferences between spectroscopic and photometric parameters and 
should be taken with caution. 

We note that star B is enriched in Al and Ba compared to 
star A, as we shown in Fig. 2. It would be tempting to assume 
that they were contributed from an unknown AGB companion 
around star B. However, it is difficult to precisely estimate the 
origin of these elements. Ba is mostly considered an s-process 
element, that is, produced by a slow neutron-capture reaction of 
relatively heavy nuclei and carried to the ISM by the winds of 
AGB stars. The isotope 7°Al could also be produced by AGB 
stars, but other possible sources are winds of massive and very 
massive stars and supernova explosions (see e.g. Martinet et al. 
2022, and references therein). In principle, a common origin for 
the two elements therefore cannot be entirely discarded. How- 
ever, this hypothetical scenario should be taken with caution 
because no AGB companion is currently detected orbiting star 
B. 


Appendix B: Scenario of the evolutionary state 


We explore in this section whether evolutionary effects might 
explain the observed chemical differences between stars A and B. 
RGB stars burn hydrogen in a shell around a inert helium core 
(Iben 1968), while red clump (RC) stars are in the stage of 
core-helium burning (Cassisi & Salaris 1997; Girardi et al. 1998). 
RGB and RC stars overlap significantly in the Teff — log g dia- 
gram (e.g. Girardi 2016) and it is hard to distinguish the two types 
of stars. Figure B.1 shows the log T,¢¢ — log g diagram for stars A 
and B (blue and red, respectively) and three (solar metallicity) evo- 
lutionary tracks corresponding to 2 Mo, 3 Mo, and 4 Mo, calcu- 
lated using the MESA isochrones and stellar tracks (MIST, Dotter 
2016)* version 1.2. In Fig. B.1, stars A and B are located between 
the tracks of 2 Mo and 4 Mo. However, it is difficult to determine 
a precise value of the mass (and evolutionary state), considering 
the mentioned overlap of RGB and RC stars in the Teff — log g 
diagram (e.g. Girardi 2016). 

Physical parameters such as mass, age, and evolutionary 
state were estimated using the latest version of the isochrones” 
package (Morton 2015). The package creates a model for 
each star and interpolates it in a grid of MIST tracks, indi- 
cating its evolutionary state through equivalent evolutionary 
points (EEPs), as described in Dotter (2016). In this scheme, 
the terminal-age main sequence (TAMS), the tip of the RGB 
(RGBTip), the zero-age core-helium burning (ZACHeB), and 
the terminal-age core-helium burning (TACHeB) correspond to 
EEP values of 454, 605, 631, and 707, respectively. We used 
as input for isochrones the high-precision stellar parame- 
ters (Terr, log g, and [Fe/H]), together with the V magnitudes 
(9.71+0.03 and 9.79+0.02 for stars A and B) and the redden- 
ing Ay. We note that the extinction in the direction of the stars is 
significant: We obtained Ay ~ 0.57 (Gaia Collaboration 2018) or 
Ay ~ 0.73 (Schlegel et al. 1998). However, the physical parame- 
ters we derived (mass, age, and evolutionary state) remain almost 
unchanged when the high-precision spectroscopic values alone 
are used (Tez, log g, and [Fe/H]). We present in Table B.1 the 
mass and age of stars A and B, estimated using isochrones, 
considering the stars in the RGB or RC phases. Depending on 
whether the stars belong to the RGB or RC phases, their phys- 
ical parameters are different and should be taken with caution. 
In Fig. B.2 we present the log Teff — log g diagram for stars A 
and B and the evolutionary tracks calculated with MIST for their 


4 https://waps.cfa.harvard.edu/MIST/ 
5 https: //github.com/timothydmorton/isochrones 
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Fig. B.1. log Teșf — log g diagram for stars A and B (blue and red, 
respectively), and three (solar metallicity) evolutionary tracks corre- 
sponding to 2 Mo, 3 Mo, and 4 Mo calculated using MIST. 


Table B.1. Physical parameters estimated with the isochrones 
package. 


Star RGB RGB RC RC 
Mass Age Mass Age 
[Mo] [Myr] [Mo] [Myr] 
+0.36 +118 +0.29 +184 
EEE 
0.27 -73 "£0.24 -133 


corresponding masses and metallicities. Blue and red correspond 
to stars A and B. The plots also show two additional tracks for 
2 Mo and 4Mo (black lines) for reference. The left and right 
panels correspond to the RGB and RC cases, respectively. 

The individual ages shown in Table B.1 strongly suggest 
that both stars belong to the RGB or that both stars belong to 
the RC, and are not a mixed pair RGB+RC. Although a num- 
ber of RGB and RC stars overlap in the log Teff — log g dia- 
gram, Bovy et al. (2014) suggested a method for distinguishing 
them. This method requires an accuracy better than 100 K and 
0.1 dex in T,¢¢ and log g. We present in Fig. B.3 the position 
of stars A and B (red and blue) in the log g - Tepp diagram and 
the cuts suggested by Bovy et al. (2014) (dashed lines) to sep- 
arate the RGB and RC regions, which is analogous to Fig. 1 
of their work. Both stars would belong to the RC region in 
this diagram. However, we caution that asteroseismic parame- 
ters for both stars could help us to distinguish their evolutionary 
stage better (e.g. Montalbán et al. 2010; Bedding et al. 2011). 
When we assume that both stars belong to the RGB or RC 
phase, we can adopt a common (average) age of 337+175 Myr or 
514+255 Myr, respectively (while a mixed pair RGB+RC would 
present an intermediate age). The common age for stars A and B 
presents a considerable uncertainty. In any case, we show below 
that chemical evolutionary effects do not explain the observed 
differences between stars A and B (for iron and other metals) for 
any combination of the considered phases. 

During the evolution of giant stars, their extensive envelopes 
mix internally processed material by fusion with unprocessed 


surface material (e.g. Salaris & Cassisi 2017). This produces a 
strong reduction of the superficial C/N ratio and lithium abun- 
dances after the FDU in the RGB phase, for example. We present 
in Fig. B.4 the C/N ratio, the lithium abundance A(Li)®, and the 
metallicity as a function of age, estimated with MIST evolution- 
ary tracks. The green bars show the Li content expected in star 
A after ingestion of planetary material (see next sections). The 
evolutionary models of stars A and B (blue and red lines) show 
a significant step or jump in C/N and A(Li), which corresponds 
to the strong reduction expected after the FDU. This step is also 
present in the evolution of the metallicity, although the jump is 
considerably lower. We note that FDU of star A occurs earlier 
than FDU of star B (the blue step occurs at a younger age than 
the red step), which corresponds to the slightly higher mass (and 
faster evolution) of star A compared to B. The low values of 
C/N and A(Li) shown by stars A and B strongly suggest that 
both stars have already gone through the FDU phase, (see, e.g., 
the right panel of Fig. B.4 for the case of A and B taken as RC 
stars). However, for the RGB case of star B (left panel, red), the 
evolutionary model predicts that the FDU will take place in the 
next few million years, while the low C/N and A(Li) values show 
that the FDU already occurred in star B. This disagreement by a 
few million years shows that the true age of both stars (for the 
RGB case) might be slightly older than estimated, or, alterna- 
tively, that star B belongs to the RC phase rather than the RGB 
phase, as suggested by Fig. B.3. 

The middle panels of Fig. B.4 show that the lithium abun- 
dances in stars A and B are lower than the expected values 
from MIST models. However, it is difficult to reach a per- 
fect agreement between observations and models because many 
effects might modify the lithium abundances. For instance, the 
modeled Li content strongly depends on the assumed value for 
V/Ucrit, Which is the initial rotation rate (compared to the critical 
one) used in the MIST models. We adopted v/v¢,;:=0.0; but if 
we had adopted v/v¢,;,=0.4, the predicted lithium abundances 
would be several orders of magnitude lower than observed. 
Charbonnel et al. (2020) showed that the evolution of Li strongly 
depends on whether rotation is included in the models (see, e.g., 
their Fig. 11). The evolutionary tracks of MIST also assume 
a starting value of A(Li) near ~3.3dex (Lodders & Fegley 
1998), while other authors adopted a lower value (3.11 dex, 
Charbonnel et al. 2020). The number of still uncertain processes 
that might modify the lithium abundances led some authors (e.g. 
Privitera et al. 2016) to suggest the use of an increased rotation 
instead of lithium, in order to detect a possible planet ingestion 
(see text above), for example. 

The evolutionary effects for the case of iron were estimated 
using MIST tracks (Fig. B.4, lower panels). The highest variation 
in [Fe/H] predicted by models during the evolution of stars A 
and B is ~0.04 dex (for the RGB or RC cases). This value corre- 
sponds to the increase in metallicity after the FDU, the step in the 
abundances shown in the Fig. B.4. However, during the evolution 
after the FDU, the predicted variation in metallicity is almost 
null (lower than 0.01 dex for both stars). Thus, we can consider 
two possibilities: First, if both stars underwent the FDU before 
(as strongly suggested by their C/N and Li abundances), evolu- 
tionary effects (<0.01 dex) clearly cannot explain the observed 
difference of ~0.08 dex. Second, if we assume that only one star 
passed through the FDU (which is highly unlikely), evolution- 
ary effects (<~0.04 dex) are not enough to explain the observed 
difference of ~0.08 dex. In this way, evolutionary effects can- 
not account for the observed differences between stars A and B, 


® We define A(Li) = log(Li/H) + 12, similar to Carlberg et al. (2012). 
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Fig. B.2. log Tes — log g diagram for stars A and B, and the evolutionary tracks calculated with MIST for their corresponding masses and 
metallicities. Blue and red correspond to stars A and B. The plot also shows two additional tracks for 2 Mo and 4M, (black lines) for reference. 
The left and right panels correspond to the RGB and RC cases, respectively. 


even assuming that only one of the stars underwent the FDU 
(which is highly unlikely). This is also valid considering that 
both stars belong to RGB, both belong to RC, or are a com- 
bination of RGB+RC. Other abundances such as [Si/H], [S/H], 
[Ca/H], and [Ti/H] were also modeled by MIST and show a step 
very similar to the steps of [Fe/H], mainly produced by a slight 
reduction of the superficial hydrogen after the FDU (Dotter et al. 
2017). 

We also performed additional experiments by calculating 
hypothetical tracks of star A. We assumed that star A was 
born with the same metallicity as star B, but considered differ- 
ent values for its mass: the same masses previously estimated 
(Table B.1) for the RGB and RC cases (dotted blue lines in the 
Fig. B.4), and also a mass of 3 Mo (dotted magenta line in the 
lower panel of Fig. B.4). An increase or decrease in the mass 
of the stars basically corresponds to a decrease or increase in 
the age when the FDU occurs (e.g., compare the dotted blue and 
magenta lines in the lower panel of Fig. B.4). In other words, 
a modified mass of the stars would basically displace the iron 
lines mainly along the x-axis in Fig. B.4 (on the other hand, a 
modified starting [Fe/H] would displace the lines mainly along 
the y-axis). Therefore, evolutionary effects does not seem to be 
enough to explain a [Fe/H] difference of ~0.08 dex between stars 
A and B (both RGB, both RC, or a combination of RGB+RC). 
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Fig. B.3. Position of stars A and B (red and blue) in the log g - Te y+ 
diagram and the cuts suggested by Bovy et al. (2014) (dashed lines) to 
separate the RGB and RC regions. The regions correspond to stars in 
the range -0.2 < [Fe/H] < 0.2. 
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Fig. B.4. C/N, A(Li), and metallicity as a function of age (upper, middle, and lower panels), estimated with MIST evolutionary tracks (blue and 
red lines for stars A and B, respectively). The green bars show the Li content expected in star A after ingestion of planetary material. The values 
measured for stars A and B are shown in blue and red (including the error bars on the y-axis). The left panels correspond to the RGB phase, and 


the right panels correspond to the RC phase. 


Appendix C: First-order estimation of the ingested 
mass 


In this section, we estimate the mass ingested by star A, if any, to 
explain its higher metallicity compared to star B. We used equa- 
tion (5) of Teske et al. (2016) (which is similar to equation (1) of 
Ramirez et al. 2011), 


A[M/H] = log (Z/X)ezMez + (Z/X)eMp . 


C.1 
ZIK aM + Mp) 1) 


where A[M/H] is the difference in metallicity, (Z/X)cz is the ratio 
of the fractional abundance of metals relative to hydrogen in the 
convective zone, Mcz is the mass of the convective envelope, 
(Z/X), is the metallicity of the planet, and M, is the mass of 
the planet. This model assumes an instantaneous accretion and 
addition of the planet to the convective envelope of the star, and 
it should be considered as a first-order approximation. Adopt- 
ing from Asplund et al. (2009) the mass fractions for the Sun, 
we have (ZX = (0.0134/0.7381) = 0.0182. We estimated the 
convection zone metallicity (Z/X)cz by scaling the solar value 
with the metallicity (similar to the procedure of Ramírez et al. 
2011; Teske et al. 2016), (Z/X)cz = 0.0182 x 107®0!5 ~ 0.0176. 
In equation C.1, we considered two reference values for 
the planet metallicity, (Z/X), = 0.10 and (Z/X), = 0.36. The first 
value of (Z/X), is similar to that of Jupiter, whose metallicity 
is estimated to be between 0.04-0.12 (Ramírez et al. 2011). For 
the second value of (Z/X),, we considered an extreme case of 
metal content inside a planet (Thorngren et al. 2016), HAT-P- 
20 b includes 600 Mg of metals inside a 7.2 Mju) (~2290 Mo) 
planet. In this last case, we roughly estimate (Z/X), ~ 600/(2290- 
600) ~0.36 for a planet with an extreme metallic content. How- 
ever, we note that these massive (and metallic) planets are 
hard to explain by current planet formation models (see e.g. 
Leconte et al. 2009; Thorngren et al. 2016), requiring a migrat- 
ing planet that accumulates all of the metal available in the disk. 
Equation C.1 also depends on the convective envelope mass 
Mcz adopted, which is a function of the stellar age. Some authors 
estimated average Mcz values near 77% of the stellar mass 


(Pasquini et al. 2007a) and ~65% (Figure 1 of Girardi 2016) 
for giant stars with masses between 1.0-1.3 Mo. Other authors 
assumed that the Mcz encompasses most of the giant star, mak- 
ing fully convective objects, that is to say, Mcz near 100% (e.g. 
Basu & Hekker 2020). Thus, for an estimate using reference val- 
ues, we considered three possible average values for Mcz: 100%, 
77%, and 35% of the stellar mass. In particular, the lower value 
of Mcz (35%) would correspond to the most favorable case for 
the detection of a metallicity signature of a planetary ingestion. 

The next step was to estimate the mass of the planet 
M, necessary to explain the difference in metallicity observed 
A[M/H]~0.08 dex. We estimated M, using equation C.1, which 
depends on the assumed values for (Z/X),, Mcz, and so on. As 
explained, we considered two possible values for (Z/X), and 
three possible values for Mcz. The combination of these val- 
ues produced six different cases with six different values for M, 
(upper panel of Table C.1). The resulting values of M, (Col. 4 
of Table C.1) show their strong dependence on the planet metal- 
licity (Z/X), and on the convective envelope mass Mcz. In these 
six cases, we assumed a mass of 2.872 Mo for star A, which 
would correspond to the RC phase (upper panel of Table C.1). 
However, if star A belongs to the RGB phase, its mass would be 
3.183 Mo, resulting in six different values for M, (lower panel 
of Table C.1). 

In addition, we took the opportunity and estimated the 
expected increase in the rotational velocity after the planet 
ingestion, AV,o; = Vyor(final) — Vv, (initial), using equation (1) 
of Carlberg et al. (2009). We also estimated the expected 
increase in the lithium abundance after the planet inges- 
tion, ALi=A(Li) fing — ACLA)jnitia, using equation (2) of 
Carlberg et al. (2012). We adopted A(Li) initia = -0.18+0.08 dex 
similar to Carlberg et al. (2012), which is the average abundance 
of giant slow rotators. Both Av,,; and ALi depend on M, and 
Mcz, and they are shown in the two last columns of Table C.1. 
In order to estimate Av,,.;, we adopted a hot Jupiter-like initial 
value of the semimajor axis (a=0.02 au) for the ingested planet. 
(For comparison, a warm Jupiter-like planet with ~0.5 au would 
imply a much higher reservoir of angular momentum, resulting 
in a significantly higher Av,,.; by several orders of magnitude). 
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Table C.1. Estimate of the required planetary mass (M,), variation in rotation (Av,,,), and lithium content (ALi) after a possible ingestion event. 


Case Mcz (Z/X)p M, AV rot ALi 
[Myup] [kms] [dex] 
Star A in the RC phase: 

1 100% 0.10 135.0 37.46 +2.10 

2 100% 0.36 31.5 8.74 +1.49 

3 11% 0.10 105.0 37.84 +2.10 

4 11% 0.36 24.2 8.72 +1.49 

5 35% 0.10 47.5 37.66 +2.10 

6 35% 0.36 11.0 8.72 +1.49 

Star A in the RGB phase: 

1 100% 0.10 150.0 37.88 +2.10 

2 100% 0.36 35.0 8.84 +1.49 

3 11% 0.10 115.0 37.71 +2.10 

4 11% 0.36 27.0 8.85 +1.49 

5 35% 0.10 52.5 37.88 +2.10 

6 35% 0.36 12.2 8.80 +1.49 
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Fig. C.1. Li abundance derived for 378 G/K giant stars from Liu et al. (2014b) as a function of T.;, and log g (left and right panels). The empty 
circles and crosses correspond to Li values and upper limits, respectively, measured for the 378 giant stars. The measured values in this work for 
stars A and B are shown in blue and red, respectively, and the predicted values after a possible planetary ingestion in star A are shown in green. 


Finally, in the estimation of Av, ; we adopted a null eccen- 
tricity (e=0), as assumed by Carlberg et al. (2009) for planets 
with unknown eccentricity. In this way, Table C.1 summarizes 
the effects after a planet ingestion, including an estimate of the 
required mass to ingest M,, together with the expected Av,,, and 
ALi values after the pollution event. 

In order to study a possible planetary ingestion, we also com- 
pared the lithium content in stars A and B with that of other 
giant stars. We present in Fig. C.1 the Li abundance derived 
for 378 G/K giant stars from Liu et al. (2014b) as a function of 
Tere and log g (left and right panels). The empty circles and 
crosses correspond to Li values and upper limits, respectively, 
measured for the 378 giant stars. We restricted the stars shown 
to a metallicity range between -0.2dex and 0.1 dex, which is 
closer to the values of the stars in this work. The measured val- 
ues in this work for stars A and B are shown in blue and red, 
respectively, and the predicted values after a possible planetary 
ingestion in star A are shown in green. The green lines corre- 
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spond to possible ALi values between 1.49 dex and 2.10 dex esti- 
mated for star A (Table C.1). Li abundances of Liu et al. (2014b) 
were corrected for NLTE effects by interpolating in the data 
of Lind et al. (2009). In order to properly compare the abun- 
dances, we also corrected for NLTE following the corrections 
of Lind et al. (2009). Fig. C.1 shows a considerable dispersion 
of Li content for a fixed value of the parameters. The Li abun- 
dances of stars A and B appear to be similar to those of other 
giant stars. Liu et al. (2014b) suggested a trend in which giant 
stars with higher Tess tend to show a higher lithium content. 
The lower Li in star A compared to star B might at least par- 
tially be explained by the higher T.¢¢ of star B. On the other 
hand, a planetary ingestion event would produce a significant 
increase in the Li content in star A (the green values shown in 
Fig. C.1), being considerably higher than current A(Li) observed 
in both stars A and B. However, Fig. C.1 shows no evidence of 
an increased Li content in star A (on the contrary, star A dis- 
plays a lower Li content than star B). In other words, there is 
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no need to invoke an ingestion event to explain the Li content in 
stars A and B. 

We also considered the possibility of a brown dwarf (BD) 
ingestion by star A. We focused on the most favorable condi- 
tions for the detection of a chemical difference, that is, by tak- 
ing Mcz ~ 35% and star A in the RC phase (mass of 2.87 Mo). 
For the metallicity of the ingested object (Z/X), (in this case, 
a brown dwarf, (Z/X)gp), we also scaled the solar value by 
the metallicity, that is, (Z/X)gp = 0.0182 x 10*°? ~ 0.0363. The 
metallicities of known BDs range between -1.0 to +0.3 dex (e.g. 
Meisner et al. 2023, and references therein), and we adopted a 
metal-rich BD with [Fe/H] = +0.3 dex. In this case, equation C.1 
implies that it would be necessary to ingest more than ~250 M jy, 
of material to reach a difference of ~0.08 dex between stars A 
and B. This is beyond the mass limits of BDs. These high mass 
values would also imply a significant increase in the rotational 
velocities and lithium content, which is not detected in the binary 
system of this work. 

The question now is when the possible ingestion occurred. 
A number of events might modify the superficial abundances 
of giant stars, such as the FDU, and therefore, it is relevant 
to consider this question. However, it is difficult to give a pre- 
cise estimate. Equation C.1 used in this work to estimate the 
possible amount of ingested mass assumes an instantaneous 
engulfment, that is, a very recent event. Similarly, equation (1) 
of Carlberg et al. (2009) (which describes the angular momen- 
tum deposition of the planet) and equation (2) of Carlberg et al. 
(2012) (which estimates the increase in lithium after the planet 
ingestion) were derived without specific time or evolutionary 
constraints. However, we can consider the following. The values 
of the C/N ratio and A(Li) strongly suggest that the giant stars 
A and B already underwent the FDU. The values shown in the 
Table C.1 show that it would be easy to detect a planet ingestion 
for a lower convective envelope mass (~35%) than for a massive 
convective envelope (~100%). This suggests that a planet inges- 
tion event would be more easily detected after the FDU rather 


than before the FDU. Similarly, Aguilera-Gomez et al. (2016) 
suggested that ingestion events would be easily detected after 
the FDU and before the RGB bump. 

For a moment, we assume the following scenario: The planet 
ingestion occurred after the FDU in star A, and several mil- 
lion years after this event, the convective envelope mixed and 
erased the additional lithium that was added by the ingestion 
(this lithium depletion is not considered by equation (2) of 
Carlberg et al. 2012), perhaps bringing the lithium abundances 
to the values currently observed in star A. However, in this case, 
it would be difficult to explain why the same mixing (which eas- 
ily restored the lithium abundance to pre-ingestion values) was 
unable to restore other metals to their pre-ingestion values. In 
other words, the massive convective envelope should diminish 
the additional lithium added by the ingestion (by mixing with 
layers with a a much lower lithium content, and perhaps also 
by reaching deeper and hotter regions in which lithium could 
be destroyed). That is, to diminish the additional lithium, sig- 
nificant mixing is required. However, in this case, it would not 
be clear why the mixing (which was significant for lithium) 
was not significant for other metals as well. Even assuming 
that the mixing erased the ingested lithium but not other met- 
als (which is highly unlikely), the values of Table C.1 for the 
case of convective envelopes of ~35% suggest that we would 
require the ingestion of massive planets (>11My,,) with an 
extremely high metallic content. As mentioned previously, these 
extremely metallic planets are very difficult to explain with 
current planet formation models and would require the accre- 
tion of almost all metals available in the disk at the time of 
planet formation. Similarly, in this case, we still have to explain 
why there is no evidence of an increased rotational velocity 
in star A compared to star B, and currently, no planets are 
detected orbiting stars A or B. Formally speaking, although 
previous arguments do not completely rule out an ingestion 
event in star A, we consider that an engulfment event is highly 
unlikely. 
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